KEK-TH-1588 



(N 



Vacuum birefringence in strong magnetic fields: 
(II) Complex refractive index from the lowest Landau level 

Koichi Hattorj*| 

Institute of Physics and Applied Physics, 

Yonsei University, Seoul 120-749, Korea 
O- 

cm; n 

q ■ Kazunori ItakuralU 

(D; 

Q i Theory Center, IPNS, High energy accelerator research organization (KEK), 

^ \ 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan and 

45 Department of Particle and Nuclear Studies, 

^ ■ Graduate University for Advanced Studies (SOKENDAI), 

1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan 
(Dated: December 11, 2012) 

l>. ; Abstract 
ov 

We compute the refractive indices of a photon propagating in strong magnetic fields on the 

7— i ; 

■ basis of the analytic representation of the vacuum polarization tensor obtained in our previous 

t-h ; 

CNJ . paper. When the external magnetic field is strong enough for the fermion one-loop diagram of 

the polarization tensor to be approximated by the lowest Landau level, the propagating mode 

• i-H ■ 

^ ■ in parallel to the magnetic field is subject to modification: The refractive index deviates from 

unity and can be very large, and when the photon energy is large enough, the refractive index 
acquires an imaginary part indicating decay of a photon into a fermion-antifermion pair. We study 
dependences of the refractive index on the propagating angle and the magnetic-field strength. It is 
also emphasized that a self-consistent treatment of the equation which defines the refractive index 
is indispensable for accurate description of the refractive index. This self-consistent treatment 
physically corresponds to consistently including the effects of back reactions of the distorted Dirac 
sea in response to the incident photon. 
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I. INTRODUCTION 

In our recent paper which we call "paper I" we analytically computed the vacuum 
polarization tensor of a propagating photon in a strong magnetic field at the one-loop level 
of a 'dressed' fermion, which is diagramatically shown in Fig. [TJ While the interaction of the 
fermion with the propagating photon (wavy lines in Fig. [T]) is considered at the lowest order 
with respect to the coupling constant, the resultant diagram contains all-order tree- level 
interactions with the external field (dashed lines) through the use of the dressed fermion (a 
double line). Such diagrams are enhanced when the external magnetic field is strong enough, 
because each insertion of the external field B yields an enhancement factor eB/m 2 = B/B c 
where m is a mass of the fermion and B c = m 2 /e is the corresponding critical magnetic field. 
The vacuum polarization tensor, or equivalently, the self-energy of a photon is a fundamental 
quantity that carries information on the change of properties of a propagating photon in 
response to the external field. 

Having obtained the analytic representation of the polarization tensor, we can now inves- 
tigate vacuum birefringence and decay of a photon into a fermion- antifermion pair, both of 
which are characteristic phenomena of "nonlinear QED effects" and should become visible 
in external fields stronger than the critical magnetic field B c . In the presence of external 
magnetic fields, the vacuum polarization tensor allows for additional terms which are absent 
in the ordinary vacuum (see Eq. ([5])). These additional terms depend on polarizations of the 
photon, and consequently lead to birefringence phenomena. Furthermore, the polarization 
tensor obtained in paper I ]] is expressed as infinite summation over all the possible pairs of 
Landau levels corresponding to those of a fermion and of an antifermion in the magnetic field, 
and show quite nontrivial behavior when the photon momentum is around the thresholds of 
the decay determined by the Landau levels. In particular, the vacuum polarization tensor 
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FIG. 1. Vacuum polarization in a magnetic field. 
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FIG. 2. Threshold structure of a photon momentum ry = g 2 /(4m 2 ) as a function of B r = 
B/B c = eB/m 2 . Curves are positions of thresholds gjj = {s/m 2 + 2^ei? + -y/ m 2 + 2(^ + n)eB} 2 
with l,n£Z, determined by the Landau levels with vanishing longitudinal momenta. The shaded 
region is densely filled with the threshold lines. The horizontal line at r 2 = 1 corresponds to the 
lowest threshold (n = I = 0). 

has an imaginary part when the photon momentum is larger than the lowest threshold which 
indicates the decay of a photon into a pair of a fermion and an antifermion, both of which 
are in the lowest Landau levels with vanishing momenta along the direction of the magnetic 
field. Moreover, each term of the infinite summation of the polarization tensor corresponds 
to a pair of Landau level indices, and thus has a unique threshold determined by the Landau 
levels. In all the terms, an imaginary part appears when the photon momentum exceeds the 
threshold specific to individual terms. 

In the present paper, we compute refractive indices and dielectric constants by using 
the analytic expression for the vacuum polarization tensor. We investigate very carefully 
their behavior around the thresholds, and how the imaginary part appears when the photon 
momentum exceeds the first threshold given by the lowest Landau levels (LLLs). In fact, 
since all the thresholds have similar structures as that of the first threshold in the vacuum 
polarization tensor, we focus in the present paper on the first threshold. In other words, we 
consider the results in the LLL approximation. This physically corresponds to working in the 
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strong limit of the external magnetic field because the first threshold is independent of the 
magnitude of the magnetic field and thus gets isolated more distantly from the higher thresh- 
olds as the magnetic field becomes stronger. This is shown in Fig. [21 Thus, our calculation 
in the present paper is valid below the second lowest threshold q 2 < (m + \Jm 2 + 2eB) . 

The refractive indices n and the dielectric constants e are defined through the dispersion 
relations for the propagating photons e = n 2 = \q 2 \/u 2 , which one obtains by solving the 
Maxwell equation with the vacuum polarization tensor included Then, one finds 

that they are in general expressed by the scalar coefficients Xi — 0, 1,2) of the vacuum 
polarization tensor (only xo remains nonzero in the vanishing field limit): 

1 + Xo m 
6± 1 + Xo + X 2 sin 2 0' U 
1 + Xo + Xi ^ 



1 + Xo + Xi cos 2 e 

where 9 is the angle between the photon momentum and the magnetic field, and the per- 
pendicular _L and parallel || indices imply polarizations with respect to the direction of the 
magnetic field. However, these expressions do not immediately provide the results we want. 
This is because the scalar functions x% are m general functions of the photon momentum 
q^ and thus contain the refractive indices (or the dielectric constants) through the disper- 
sion relation. Thus, to obtain the correct values for the refractive indices (and dielectric 
constants), one has to treat the equations self-consistently. We will see that such a self- 
consistent procedure gives a large effect when the deviation of the refractive indices (or the 
dielectric constants) from unity is large. 

As is well explained in standard textbooks on optics, refractive indices have clear physical 
meaning: Real part of the refraction index provides a phase velocity i> p hase = 1/^reai of a 
propagating photon and imaginary part is called an extinction coefficient k = n imag for the 
photon propagation. They appear in a phase factor and a damping factor of the photon 
field *if(t : x) (see Appendix D in Ref. jl| for detail): 

V(t,x) <xexp{-iu(t-ir^ sc q- x)} exp{-UKq-x} , (3) 

where q is a unit vector directed to the photon propagation. Thus it is natural to define 
decay length in which the intensity of the photon field I(t, x) oc \^(t, x)\ 2 falls off in a strong 
magnetic field, 
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When the refractive index has an imaginary part, the decay length takes a finite value, which 
also indicates that the photon decays within a finite lifetime. Notice that these quantities 
depend on the photon energy u and the angle 9 between the magnetic field and the photon 
propagation direction. In fact, they are phenomenologically important when we apply our 
results to realistic situations such as in heavy-ion collisions and high-intensity laser where 
strong magnetic fields are realized only in a small space-time volume. Namely, if decay length 
and lifetime are small enough compared to the spatial extent and time duration of strong 
magnetic fields, then the decay process will be dominant. Otherwise only the birefringence 
phenomena will be seen. 

The present paper is organized as follows: in the next section, we briefly review the 
main results obtained in paper I. In particular, we give the explicit form of the vacuum 
polarization tensor in the LLL approximation which is used in the present paper. Since the 
perpendicular polarization of the dielectric constant is not modified (ej_ = 1) in the LLL 
approximations, we will discuss only the parallel polarization en in the following sections. 
In Sec. Ill II we discuss general aspects of the dielectric constant in the LLL approximation, 
which is followed by detailed description in Sec. IIVI There, we compare results from different 
methods in solving the self-consistent equation for the dielectric constant. Then, in Sec. |V] 
we show the refractive index of the parallel component. Since the refractive index is trivially 
related to the dielectric constant, the results are qualitatively the same, but still, it will be 
useful when we apply our results to phenomenology. The last section is devoted to summary 
and outlook for future applications in laser physics, heavy-ion physics, and compact stars 
with high magnetic fields. In Appendices, we provide some supplementary information such 
as the dependence of dielectric constant on the propagation angle and the magnetic field 
strength, and a brief comment on other higher Landau levels. 



II. VACUUM POLARIZATION TENSOR IN STRONG MAGNETIC FIELDS 

In this section, we briefly summarize the main results obtained in paper I We also 
provide notations necessary for the calculation in the present paper. As we mentioned in 
Introduction, we work in the strong field limit and thus we present here the explicit analytic 
1 The other contributions from the higher Landau levels barely modify the perpendicular polarization e± 

in the kinematical region below the lowest threshold. This will be briefly mentioned again in associate 

with Fig. H 
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expression of the vacuum polarization tensor in the LLL approximation. 



A. Vacuum polarization tensor and dielectric constants 

Suppose that the external magnetic field is applied along the third axis of spatial coordi- 
nates in the negative direction so that eB 3 = \e\B > for an electron (£> 3 = —B is the third 
component of the magnetic field vector B l , and e = — |e| is negative for an electron). The 
photon momentum q^ is now decomposed into longitudinal and transverse components with 
respect to the magnetic field, namely = (g°,0,0,g 3 ) and q^_ = (0, 9 1 , g 2 , 0). Correspond- 
ingly, the metric tensor t]^ = diag(l, —1,-1,-1) is also decomposed into longitudinal and 
transverse subspaces r)u U = diag(l, 0, 0, —1) and rf± = diag(0, —1, —1, 0). 

In the absence of the full Lorentz symmetry, the vacuum polarization tensor U^(q\\, q±; B T ) 
allows for two additional terms. Using the projection operators Pj^ satisfying q^P^ = 0, a 
gauge-invariant tensor structure is given by (2 6 1 (see also Refs. jl, 8] for details), 

iC(?ii, ?x; B t ) = - (xoPr + XiPT + X*Pf) , (5) 

PF = frr-M, PT = ihf - Pr = qivT-g1ql, (6) 

where we suppressed arguments of Lorentz-scalar coefficient functions Xi = 0, 1,2). All 
the three coefficient functions depend on the magnetic field B T = B/B c , but only xo survives 
nonzero in the ordinary vacuum with the vanishing magnetic field. 

Plugging the general form of the polarization tensor (jSJ) into the Maxwell equation for the 
propagating photon, we can solve it to obtain the dispersion relations for physical modes 
(see Sec. 3 and Appendices C and D in paper I |l|). As mentioned in Introduction, the 
dielectric constants e or the refractive indices n are defined by 



n 2 



\Q 



2 

(7) 



where the photon momentum is q^ = (u, q). Then, one finds two different values e± and eii 
for the dielectric constants as shown in Eqs. ([T|) and fl2]). There, we have taken a coordinate 
system in which a photon is propagating in the y = plane so that the momentum vector 
g M is represented as g M = (ou, \q\ sin(7r — 9), 0, \q\ cos(7r — 8)) with 9 being the angle between 
the direction of the external magnetic field and the momentum of a propagating photon. 

Since two dielectric constants e± and ey, and consequently, two refractive indices n± and 
tin, are in general not equal to unity and different from each other, we call this phenomenon 
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vacuum birefringence after a similar phenomenon in dielectric substances. Notice also that, 
due to the violation of the Lorentz symmetry by an external magnetic field, the dielectric 
constants explicitly depend on the zenith angle 9. Nevertheless, the system maintains boost 
invariance in the third direction, and thus photon propagations in the directions at 9 = 0, 7r 
are special. Indeed, substituting these angles into Eqs. ([T]) and (j2J), we find that both of 
them become unity e±(9 = 0,7r) = e\\{9 = 0, ir) — 1 and that there is no effect from the 
magnetic field. 

As we will discuss later and also discussed in detail in paper I the scalar coefficients 
can become complex when the photon momentum exceeds some threshold values. In that 
case, the corresponding dielectric constants and refractive indices are also complex. Since 
the scalar coefficients Xi are parts of the self-energy of the photon, emergence of an imaginary 
part indicates the decay of a photon into a fermion-antifermion pair. 



B. Scalar coefficient functions in the LLL approximation 

Let us outline how to obtain analytic representations for the scalar coefficient functions 
Xi in the external magnetic fields. In paper I, we computed Xi by Schwinger's proper 
;ime method which is a very useful technique for the calculations under external fields 



9] 



. Each coefficient function was first expressed as a double integral with respect to two 
proper times, each of which is associated with the fermion propagator composing the fermion 
loop (see Fig. [1]). Then, the integration was analytically performed after we made rapidly 
oscillating integrands tractable by expanding them with known special functions. The results 
we obtained are represented as infinite summation with respect to two indices, which are 
afterward identified with the Landau levels. In particular, the first term in the double infinite 
sum corresponds to the contribution of the lowest Landau levels, which was confirmed to 
agree with expressions obtained by other approximate methods in the strong field limit 



101 ] . We will investigate this contribution in detail in the present work. See paper I for 
details of the calculation and the explicit analytic forms of all the scalar coefficient functions 

Q. 

Here we just show the explicit form of the coefficient functions from the lowest Landau 



levels. Below, we use the following notations for dimensionless kinematical variables: r 



qy(Am ) and r\ = q\/{Am ) = — \qj_\ /(Am ). We also use a shorthand notation: r\ = 



—2r\/B r . By taking the first term in the sum with respect to the Landau levels (see 
Sec. 4.2 in paper I |l| for details), we find that the coefficient functions \% °f t ne vacuum 
polarization tensor survive only for % — 1: 



x^ LL ( 
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where a = e 2 /(4ir) and /oA( r y) * s a function responsible for the threshold behavior at the 
lowest Landau level. Depending on the value of rjj, the piecewise expression of /qa^h) i s 



given as 



i 



/rf(rf-l) 
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/rf(l-rp 
1 



hi 



arctan 



^j(i-j) 
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(rjj < 0) 

(o < rj < r 

(1 < rt) 



(10) 



Notice that this function has a singular behavior at rjj = 1 and that rjj = 1 corresponds 
to the lowest threshold for the decay of a photon into a fermion-antifermion pair. This is 
confirmed by the following observations. Firstly, the function /oA^y); an d thus the scalar 
function has an imaginary part in the region, r | > 1. Next, as briefly commented in 

Introduction and shown in paper I pj, each term in the infinite summation of the coefficient 
functions is specified by two integers n and £, and acquires an imaginary part as the photon 
momentum increases beyond the threshold value rjj = {y/l + 2£B T + ^1 + 2(£ + n)^} /4. 
This representation of the threshold suggests that £ and £ + n can be interpreted as the 
Landau levels of a fermion and an antifermion. From these, one finds that the condition 



1 is realized by the lowest Landau levels 



n 



0, and thus emergence of an imaginary 



part in rjj > 1 implies the decay of a photon into a fermion-antifermion pair in the lowest 
Landau levels. Note also that /oA( r y) does not depend on B r , reflecting the fact that the 
location of the threshold for the lowest Landau level, rjj = 1, is independent of B r . Therefore, 
within the LLL approximation, the coefficient function x\ LL linearly scales as the magnetic- 
field strength B r increases when r\ B r (see Eq. (j5J)), which will be also seen in Fig. H] 
when the photon momentum is below the lowest threshold. Finally, it should be emphasized 
that rjj = 1 is the only one singular point of the function /oA( r y)- Indeed, it is continuous 
at r« = with a finite value /qa(0) — 2, without any singularity. In Fig. [31 we show the 
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FIG. 3. The scalar function Xi ■ blue and red lines show real and imaginary parts of X\ 

behavior of the scalar function x\ LL around the threshold rj? = 1. The other thresholds from 
higher Landau levels have qualitatively the same structure as shown in Fig. El 

III. GENERALITIES ON DIELECTRIC CONSTANT IN THE LLL APPROXIMA- 
TION 

By using the analytic expression of the vacuum polarization tensor, we can compute 
the dielectric constants given in Eqs. ([1}) and fl2]). In this section, we discuss some general 
aspects of the dielectric constants and the difficulties in computing them from the coefficient 
functions Xi- A complete description of the dielectric constant within the LLL approximation 
will be presented in the next section. 

First of all, we emphasize that we have not yet specified any dispersion relation for the 
external photon momentum in the scalar coefficient functions x%- F° r an on-shell photon, the 
dispersion relations should be determined by Eqs. ([[)) and ()2]), which do contain the dielectric 
constants on the right-hand sides through the photon momenta, rjj and r\. According to 
the definition of the dielectric constant ((?]), those photon momenta are written as, r| = 
u 2 {l — e cos 2 9) and r\ = — eu 2 sm 2 9, commonly for and e\\, where we introduced a 
scaled photon energy, u 2 = u 2 /(4:m 2 ). This fact indicates that we have to solve these 
relations in a self-consistent way with respect to the dielectric constant appearing on the 
both sides. As we will find later, there is a large deviation from the massless-type dispersion 
relation indicating that r 2 + r 2 ± ^ and u ^ \q\. We will see that a self-consistent treatment 
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plays an important role in such a case. In particular, we will investigate very carefully how 
the presence of the threshold and an imaginary part in x% affects the results, which can be 
well studied in the LLL approximation. 

Within the LLL approximation, only Xi is nonzero as shown in Eqs. (jSj) and (jH]). It 
immediately follows from Eq. ([T]) that the dielectric constant in one of the polarization 
modes is intact, e± = 1. This is understood from following two observations. First, as shown 
in Appendix D in paper I an electric field induced in this polarization mode oscillates 
perpendicularly to the external field. Second, since the discretized transverse momentum of 
a fermion is fixed to the lowest Landau level, motion of a vacuum pair excitation is restricted 
only along the external magnetic field. Because the direction of the possible dipole excitation 
is perpendicular to the electric field accompanying this polarization mode, this mode does 
not induce an electric dipole excitation. Therefore, its propagation is not modified. On 
the other hand, an electric field accompanying the other polarization mode has a parallel 
component oscillating along the external magnetic field. Thus, an electric dipole can be 
induced as a response of the Dirac sea, which results in a modification of the dielectric 
constant en shown in Eq. (j2J). 

Inserting Xo — X2 — into the general expressions ([I]) and (j2J), we find the dielectric 
constants in the LLL approximation: 



6^ = 1, (11) 

u.,,- :f , „. 1 + Xi^(j,rl;i? r ) 

e,, (u, 0, B T ) - - — t t t / 2 „2 . p \ „ a ■ ( 12 ) 



Here, we have explicitly shown the arguments of ej^ 1 " and Xi LL - The dielectric constant e^ LL 
depends on a scaled photon energy u 2 = uj 2 / (4m 2 ) and angle 9 separately, because of an 
explicit ^-dependence on the right-hand side. It should be noticed that photon's squared 
momenta on the right-hand side can be rewritten in terms of the scaled energy Cj so that 
they contain the dielectric constant ejj" 1 " 1 " according to the definition ([7]): 



2 -Cu 2 (l-e\ LL cos 2 6), (13) 
I = -Co 2 e\ LL sm 2 e. (14) 



The expression of the dielectric constant (1121) has to be solved with respect to en. We 
will discuss in detail that this procedure plays an important role when modification of the 
dielectric constant becomes sizable. 
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Strong field limit (LLL) 

Soft photon & weak field limits 
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FIG. 4. Independence of the scalar functions Xi (* = 0, 1) 2): three black lines show results from 



numerical integration 
limits, respectively. 



Ill ], and red and blue lines show the weak-field 



5[ and strong-field (LLL) 



The vacuum birefringence in the LLL approximation bears an analogy to 'uniaxial' bire- 
fringent material, of which optical axis corresponds to a preferred orientation provided by the 
external magnetic field. In this analogy, the polarization mode having the dielectric constant 
en could be called " extraordinary mode", while the other one having e± " ordinary mode". 
However, this analogy holds only in the LLL approximation because the perpendicular part 
e_i_ in general depends on the angle 9. Thus, if one looks at the wide range of kinematical re- 
gion, the vacuum birefringence is much more complicated than the birefringence in uniaxial 
materials. 



A. Below threshold: Accuracy check of the LLL approximation 

Consider first the region below the threshold rjj < 1. As long as the photon momentum 
is in this region, one can numerically perform the double integration in a safe way! 2 ! Since 
the original double integrals are directly evaluated (without being expanded with respect 
to special functions), the numerical results should contain all the contributions from Lan- 

n 

dau levels in our analytic expression obtained in paper I [1]. Such numerical integration 
2 Each function x% contains a double integral with respect to two proper-time variables T\ and T2, and the 
integration over r = eB(j\ + 72)/2 in general suffers from strong oscillation of the integrand. However, 
it can be shown that for r\ < 1 an integral contour in the complex r-plane can be safely moved to the 

nn 

imaginary axis to obtain a well-behaved rcal-valujj integral which allows for numerical evaluation [3, LL1[ . 



was performed in Ref. (llj some years ago. Thus, we can check the validity of the LLL 
approximation by comparing the analytic results fl9]) & f fTOj) with the numerical results. 

In Fig. HI we compare the LLL results & (TTU]) for the coefficient function Xi LL with 
Xi from the numerical integration below the threshold. The results of numerical integration 
(three black lines) shown in Fig.H]were reproduced by ourselves, and agree with the previous 



results in Ref. up to small invisible numerical errors. We also show an approximate 
expression for soft photon and weak-magnetic- field limit j^j. The photon momenta are 
taken to be small = — r\ = 0.01 below the threshold, as an illustration. In fact, this 
is in the 'soft gamma-ray' regimqj which corresponds to typical burst spectrum from the 



soft-gamma-ray repeaters as known as magnetars. The focus of Ref. (ll| was the radiation 
from magnetars, which is the reason why they discussed the region below the threshold. The 
figure shows that the result of xi from numerical integration agrees with analytic expressions 
in the weak and strong field limits. In particular, the LLL approximation in xi works well 
for B r > 10, and the values of xo and X2 are very small compared to unity, both of which 
confirm the use of the LLL approximation (IE]) ^ (EJ) in strong magnetic fields. As mentioned 
below Eq. ( ITUl) . we find that xi linearly scales as B r increases in the strong magnetic field 
limit. 

As long as the photon momentum is below the threshold rjj < 1, it is not technically 
difficult to obtain the dielectric constant and the refractive index. They have only real 
parts, and we can use the coefficient functions obtained by numerical integration. This 



was already performed in Ref. |llj . where the results show that the refractive index in the 
parallel mode, n\\, deviates from unity as the photon momentum approaches the threshold 
rjj = 1 from below. (Essentially the same result is shown in Fig. [5] for the dielectric constant 
in the parallel mode ey.) Note that this method is valid for arbitrary strength of the external 
magnetic field, whose validity is however restricted to the kinematical region rjj < 1. It was 
also shown that deviation of n\\ from unity becomes larger as the magnetic field strength 
increases beyond the critical field strength B T = B/B c > 1. On the other hand, the refractive 
index is barely modified in the weak field region below the the critical field strength B T < 1, 
when the photon momentum is below the threshold rjj < 1. We will compare the results 
from the numerical integration with those from the LLL approximation in later sections. 
The refractive index in the perpendicular mode, n±, was also examined for r?, < 1 in 



3 Electron mass is assumed here for the fermion mass "m" in r?, and r\. 
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Ref. 



11] , and the result was reproduced by the present authors. The results do not show 



sizable deviation from unity, even when the magnetic field strength becomes larger than 
the critical field strength and when the photon momentum approaches the lowest threshold. 
Figure H] shows that Xo an d X2 are smaller than unity by two orders in magnitude, and thus 
they bring in only tiny modification in the dielectric constant in Eq. flTJ, and also in the 
refractive index n±. These results are quantitatively consistent with the dielectric constant 
(II ip in the LLL approximation. 



B. Beyond threshold: Appearance of imaginary part 

Before we present complete self-consistent treatment of the dielectric constant, we explain 
how the dielectric constant behaves around the threshold in a very rough estimate. First of 
all, let us briefly recall the behavior of Xi LL given in Eqs. (j9j and ( ITU]) around the threshold. 
As described below Eq. ([TO]) , the function Xi LL has a singularity at the threshold rjj = 1. 
When rjj goes beyond 1, it starts to have an imaginary part. The real and imaginary parts of 
x\ LL diverge when the photon momentum approaches the threshold r 2 = 1 from below and 
above, respectively (see Fig. [3]). Whereas the emergence of an imaginary part in Xi LL directly 
affects the dielectric constant e^ hh (namely, e^ Lh starts to have an imaginary part above the 
threshold), we should notice that the dielectric constant takes a real finite valuej at the 
threshold owing to cancellation of the divergences between the numerator and denominator 
in Eq. ( 1T2|) . In fact, we find a simple expression for the value of en at the threshold (for 
cos 9 ^ 0): 

lim eiifrj?) = -—_=*. , (15) 

rf-H±0 " 11 COS 2 6 " m ' V 7 

when r| approaches the threshold from both below and above. Note that the limiting value 
( !T5|) is valid even if we include contributions from other levels [I 7^ 0, n / 0). This is 
because the divergent contribution to xi & t the threshold dominates any finite contribution 
from other levels (That is why we do not put "LLL" in the above equation). The same is 
true for the other higher thresholds: the limiting value of en at each threshold comes from 
the contribution in x% which gives divergence at that point. 

4 This is trivially true as long as cos 2 9 is not zero, 9 7^ ±7r/2. In two particular cases at 9 = ±7r/2, the 
dielectric constant is given by ey = l + Xi LL with r 2 = w 2 , which directly reflects the behavior of Xi LL , and 
is also real below the threshold. These angle- dependent behaviors will be again mentioned in Sec. IV Al 
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eiim 2 
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FIG. 5. Dielectric constant in the LLL approximation: blue and red lines show the real and 
imaginary parts of the right-hand side of Eq. (I12p with eJj' LL = 1 in Xi LL (called solution (I) in 
Sec. IIV A|) Filled backgrounds indicate the kinematical regimes below and above r 2 = 1 whose 
boundary is given by Co 2 = 1/(1 — cos 2 6) assuming efj" 1 " 1 " = 1 in r|. Green dots show the numerical 
computation of the real part in r| < 1. 

Figure shows a plot of the dielectric constant f fT2"j) obtained with xi shown in Eqs. 
and (fl"0|) as a function of the photon energy squared Cj 2 . We have taken 9 = 7r/4 just for 
illustration (then, eii m = 2), and a relatively large value for the magnetic field B r = B/B c = 
500 so that the LLL approximation is appropriate. The dielectric constant which implicitly 
appears on the right-hand side of Eq. ( ji~2l) is, at this moment, taken to be one. Then, 



inserting e^ hh = 1 into r 2 given by Eq. ( fl3l) . one finds that the threshold condition r| = 1 
is now represented in the energy space aqj Co 2 = 1/(1 — cos 2 9) = 2 for 9 = 7r/4. The region 
uj 2 < 2 filled with light blue color in Fig. |5] corresponds to the one below the threshold, while 
the region u 2 > 2 with pale red, beyond the threshold. 

First let us see the region below the threshold. Green dots correspond to the dielectric 
constant ei| shown in Eq. (T5]) obtained on the basis of the numerical integration of Xo and 
Xi performed below the threshold [] (the same result as obtained in Ref. 11]). On the other 
hand, a blue line corresponds to the result from the LLL approximation. As we already 

5 Note that this equation is true only in the naive prescription where we assume = 1 in r| . 

6 To remove the divergence of xo ; we have employed a cut-off and renormalization condition described in 



the Appendix B in paper 



.a 
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discussed in the previous subsection, both results grow with increasing photon energies. We 
will discuss the deviation between these two in the next subsection. The imaginary part 
shown with a red line is exactly zero in this region. 

Next, look at the region beyond the threshold, u 2 > 2, in Fig. It is clearly seen that the 
dielectric constant has a sizeable imaginary part. This contribution is never described by the 
result from numerical integration which is valid only below the threshold. We also find that 
the real part shows nontrivial behavior: starting from Refejj" 1 " 1 "] = 2 at u 2 = 2, it decreases, 
and even becomes smaller than unity. Considering the deviation between the results from 
the numerical integration and the LLL approximation in the region below the threshold, we 
expect these results (without self-consistent treatment) will not be quantitatively accurate. 
We will indeed see in the next section that the solution will be greatly modified in the self- 
consistent treatment. What is sure at this moment is the presence of the imaginary part in 
this region. 

C. Necessity of self-consistent description 

Let us come back to the results below the threshold. We have seen that there is deviation 
between the results from the numerical integrations (green dots) and the LLL approximation 
(a blue line), and that it becomes larger as the photon energy approaches the threshold. 
Since the results grow with increasing photon energies, one can also say that two results 
show a good agreement when the dielectric constant is close to unity, but deviate from each 
other when it becomes large. Clearly, the analytic result with a large value of en > 1 is not 
consistent with the prescription assuming en = 1 on the right-hand side of Eq. f)12p . On the 
other hand, to obtain the result from numerical integration, Eq. ([2]) was solved with respect 
to u with its dependence on the right-hand side taken into account. Considering the facts 
that the result from numerical integration should be accurate in this region and that the 
magnetic field is strong enough for the coefficient functions to be approximated by the LLLs, 
we can easily conclude that the deviation of two results nearby the threshold originates from 
the "inconsistent" treatment of the dielectric constant in the naive prescription. 

There is another evidence for the importance of the self-consistent treatment. Note 
that the analytic expression of the limiting behavior shown in Eq. (I15p is obtained from 
the divergent behavior of Xi when the longitudinal momentum approaches the threshold, 
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r 2 — >■ 1. Inserting the limiting value e lim into r 2 = u) 2 (l — ey cos 2 6 1 ), we however notice that, 
as long as Q is finite, r 2 approaches zero when the dielectric constant approaches the limiting 
value e.. , 

lim r 2 (en) = . (16) 

Therefore, Eqs. ( fT5l) and ( 1T61) cannot be consistent with each other. This occurs when the 
dielectric constant is large, and must be resolved if we treat it self-consistently. 

Then, what kind of physics is involved in the self-consistent description of Eq. (j!2p ? It is 
helpful to recall a microscopic picture of the propagation of a photon in an ordinary medium. 
Variation of the dielectric constant or the refractive index is induced by the response of 
the medium to the incident photon. In the present case, medium is the vacuum which is 
filled with fermions in the Dirac sea, and the incident photon creates a polarization in it 
and, when its energy is large enough, even an on-shell fermion and antifermion pair which 
might be, alternative to the photon decay, seen as 'photoelectric effect' by absorption of 
the incident photon. These microscopic processes will generate 'back reactions' leading to 
screening and damping of an incident photon field. Notice that these effects are embedded 
in the equation as modification of the external photon momenta r 2 and r\ through the 
implicit dependences on the dielectric constant (see Eqs. (fl3l) and (1141) ). Therefore, what is 
necessary for examining these effects is, in technical terms, to solve Eq. ( fl2l) self-consistently 
with respect to the dielectric constant appearing on the both sides. 

IV. SELF-CONSISTENT DESCRIPTION OF DIELECTRIC CONSTANT IN THE 
LLL APPROXIMATION 

In the previous section, we have observed that there is deviation between the results 
from the LLL approximation (with a naive prescription) and the numerical integration as the 
photon energy approaches the threshold r| = 1 from below, and that there is a contradiction 
between the limiting value of the dielectric constant at the threshold (1151) and the behavior 
of the squared photon momentum (fl6i) . As already suggested in Sec. IHI CI it will turn out 
that these two problems originate from "inconsistent" treatment of Eq. (Tl2l) with respect to 
the dielectric constant ^ hh > an< ^ ^ na ^ they are no ^ P resen t i n fully self-consistent description 
of Eq. (|12j) . In this section, in order to be convinced of the importance of fully self-consistent 
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description, we increase step by step the levels of accuracy in treating Eq. ( fl2|) . eventually 
leading to the completely self-consistent solution. 

A. Three steps towards fully self-consistent description 

As discussed in Sec. IIII C\ Eq. ( I12p contains nontrivial effects such as screening and 
damping of an incident photon field, and they are correctly described in a self-consistent 
description. In order to entangle these effects and understand physics consequences arising 
from them, we solve Eq. ( fl2|) at three different levels of accuracy (I), (II) and (III) in treating 
the coefficient Xi LL - We also compare the results of these three methods with the result 
from numerical integration (Num.) which is valid only below the threshold. 

(I) "inconsistent" solution: Instead of solving Eq. ( fl2l) with respect to ej[ LL , we just 
assume ejj" 1 " 1 " = 1 in the coefficient function Xi LL 011 t ne right-hand side. This should 
be a good approximation when the deviation of e^ Lh from unity is small enough. The 
result of this treatment was already shown in Fig. [5j 

(II) partially self-consistent solution: We allow for a real part of e| LL on the right- 
hand side of Eq. ( I12p . Namely, the equation for the real part Re[e| LL ] is solved self- 
consistently, but the imaginary part Im[e| LL ] is just given with the scalar function 
Xx 1X1 which ejj" 1 " 1 " is replaced by the solution of the real part. Thus, this treatment 
gives a partially self-consistent solution and should be a good approximation when 
the imaginary part is small enough. In particular, this treatment is valid below the 
threshold where e| LL does not have an imaginary part. 

(III) fully self-consistent solution: We treat e^ LL in x\ LL as a complex number, and 
solve the real and imaginary parts of Eq. ( fl2l) fully self-consistently to find the de- 
pendence of eY" L on a photon energy. This is the most accurate description of the 
dielectric constant in the LLL approximation. 

(Num.) self-consistent numerical solution below the threshold: As already 
mentioned in the previous sections, one can directly evaluate the scalar functions 
Xi {i = 0,1,2) by numerical integration as long as the photon momentum is below 
the threshold r?, < 1. One can also solve Eqs. ([T]) and ([2]) self-consistently. However, 



17 



since this is possible only below the threshold, the dielectric constants are real. This 
treatment should be compared with (II) and (III) below the threshold. 



B. Threshold condition in self-consistent description 

Before we present results for the dielectric constant in the self-consistent description, let 
us revisit the threshold condition rjj = 1. We discussed that the threshold condition rjj = 1 
gives the threshold photon energy u 2 h = 1/(1 — cos 2 #) (see Fig. [5]), if en = 1 is assumed in 
rjj. In this case, the region below the threshold is simply represented by Co 2 < u) 2 h as shown 
in Fig. [5] (where u 2 h = 2 for 9 = it /4). However, it should be noticed first that, if we keep 
the real part of eii in rjj, the threshold photon energy is given in an en-dependent form as, 



and this condition is inversely solved as, 



The above relation ffl8|) separates the kinematical regions below and above the threshold 
of photon decay. Therefore, in the ey-cu 2 plane, the threshold line is not just a vertical 
line (as shown in Fig. E]), but is a curve given by Eq. ( TT8"j) . Inserting the limiting value 
e lim (EES]) into Eq. (IPT]) . we find that the photon energy approaches infinity, Cj — > oo, as the 
dielectric constant approaches the limiting value. Equivalently, Eq. ([18]) shows that the 
dielectric constant does not reach the limiting value e lim at r 2 (ey) = 1, as long as the photon 
energy Cj is finite. This observation implies that the dielectric constant could stay in a real 
value when Cj increases from zero to infinity, and thus that the energy dependence of the 
dielectric constant would be significantly modified from the one in Fig. obtained without 
any self-consistent treatment. 

C. Back reactions in complex dielectric constant 

Now let us finally show results from the self-consistent description. Figure is a compi- 
lation of all the results of the real and imaginary parts of the dielectric constant obtained 
with four different treatments from (I) to (III) and (Num.). Upper and lower panels show 
the real and imaginary parts, respectively. We have taken the same values for B T and 9 as 




(18) 
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Stable branch (I) Stable branches (II & III) • Stable branch (Extended Num.) 

— Unstable branch(I) :::::::: Untable branches (II) Unstable branch(III) 



FIG. 6. Dielectric constant for 9 = tt/4 and B r = 500: a shift of the threshold in solution (III) 
from that in solution (I) is indicated with a green stripe in the lower panel. 



in Fig. [5j namely B T = 500 and 8 = ir/A. The threshold boundaries specified by Eq. (|T8|) 
and r^(e\\ = 1) = 1 are indicated with green solid and dashed lines, left (right) of which is 
filled with light blue (pale red) and corresponds to the region below (above) the threshold. 
Long dashed lines (blue and red) are the "inconsistent" solutions from (I), short dashed lines 
(purple and red) are the partially self-consistent solutions from (II), and finally thick bold 
lines (blue and red) are the fully self-consistent solutions from (III). These are compared 
with green dots which are the results of numerical integration (Num.). Notice that this is 
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the same result as shown in Fig. [5j but is now extended to larger Co 2 because the threshold 
line is modified from a vertical line to a curve and the region below the threshold (filled with 
light blue color) extends to infinite Co 2 . 

1. Screening of incident light 

Let us first consider the change from the "inconsistent" solution (I) shown as blue and red 
long-dashed lines, to the partially self-consistent solution (II) shown as a blue solid line and 
purple and red short-dashed lines. Look at the blue solid line in the upper panel, which is a 
common solution in treatments (II) and (III). Also, it coincides with the green dots which 
are from the numerical integration (Num.). Thus, as expected, the treatment (II) already 
gives a final result in this region below the threshold. Since this line is not accompanied by 
an imaginary part at any photon energy, we call this stable branch. In particular, we now 
understand that the origin of the deviation found in Fig. [5] between the results from the 
LLL approximation with the naive prescription (blue long-dashed line) and the numerical 
integration (green dots) is the incomplete treatment of the equation (112j) . and that the 
value of the dielectric constant in the naive prescription is reduced (i.e., screened) due to 
back reactions. In addition to the stable branch, we found other two solutions (purple and 
red short-dashed lines) in the region, T 2 (e||) > 1. Namely, in this partially self-consistent 
treatment (II), we obtained the dielectric constant as a three- valued function of the photon 
energy Co. The latter two branches in the region, Tjj(e||) > 1, smoothly connects to each 
other. In contrast to the solution in the region, Tjj(e||) < 1, these branches are accompanied 
by imaginary parts as shown in the lower panel, and thus we call them unstable branches. 

As described above, the green solid line divides the kinematical regions f 2 (e||) < 1 and 
r|(e||) > 1, while the green dashed line ?" 2 (e|| = 1) < 1 and rjj(e|| = 1) > 1. Noting a 
correspondence between these lines, we find that the partially self-consistent solution (II) 
still maintains a topological structure of (I) as correspondences between stable branches 
and unstable branches, respectively: the blue long-dashed (red long-dashed) line is modified 
to the blue solid (purple and red short-dashed) line in the solution (II). Location of the 
connecting point between stable and unstable branches, where the dielectric constant takes 
the limiting value e lim , has moved away from Co 2 = 1/(1 — cos 2 6) in (I) to Co 2 = oo in (II). 

The lower panel in Fig. [6] shows an imaginary part of the dielectric constant associated 
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with the unstable branches shown in the upper panel. Since the partially self-consistent 
solution (II) gives the double-valued unstable branch, an associated imaginary part is also 
a double-valued function of the photon energy. As shown with a green stripe in the back- 
ground, the threshold of the photon decay shifts upward because of the modification of the 
real part incorporated in r|(e||). 

2. Damping of incident light 

Now let us see the change of solutions from (II) to (III). In partially self-consistent 
treatment (II), we treated the dielectric constant in Xi LL as a rea l value. However, since 
the unstable branch has a large imaginary part (see purple short-dashed line in Fig. [6]), we 
should not assume that the dielectric constants in r?, and r\ are real values. Physically, 
this indicates that we have to self-consistently take into account the damping of an incident 
photon due to decay into a fermion-antifermion pair in the external magnetic field when the 
decay rate is sizable. 

We thus solve simultaneously the equations for real and imaginary parts of Eq. (fl2"|) 
without any assumption. This is the result of treatment (III). Maintaining both the real 
and imaginary parts of the dielectric constant in the photon momenta and r\, we find 
that the structure of the branch in the kinematical region above the threshold is significantly 
modified. The most remarkable modification of the structure is that an unstable branch in 
the real part, indicated with a red solid line, directly connects to the stable branch below 
the threshold. According to this, the dielectric constant reduces from the three-valued 
function to a double-valued function of the photon energy. This edge of the unstable branch 
close to the connecting point is still accompanied by an imaginary part as shown in the 
lower panel, and thus is a state unstable to the decay. We find that a magnitude of the 
imaginary part is diminished in (III) compared to those of the imaginary parts obtained 
in the other two cases, (I) and (II). This result can be intuitively understood as follows: 
the magnitude of the imaginary part should be suppressed in the self-consistent treatment, 
since the amplitude of the incident photon damps due to the decay. Indeed, the amplitude 
of the photon propagation is given by the complex refractive index related to the dielectric 
constant as shown in Eqs. fl3]). In Sec. |V] we will also show a complex refractive index 
obtained from the dielectric constant shown in this section. 
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FIG. 7. Dielectric constant for 9 = ir/4 and very strong magnetic field B r = 5000: each line shows 
the corresponding quantities in the legend of Fig. [6j 

D. Results in stronger magnetic field 

We shall see a dielectric constant in a stronger magnetic field. Figure [7] shows the dielectric 
constant at much larger value of B v = 5000. Each line shows the corresponding quantities 
indicated in the legend of Fig. |6j Compared to the results for B r = 500 in Fig. |6l one finds 
that the modification of the real part in (I) is larger, in particular in u 2 > 2, with a nontrivial 
wavy behavior. However, as shown in red solid line of treatment (III), it is suppressed if 
we self-consistently incorporate a screening effect on the incident photon field caused by the 
induced vacuum polarization. Consequently, the real part eventually behaves similarly as 
in the previous case with B T = 500, but its magnitude shows moderate enhancement. A 
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significant effect of the self-consistent treatment is also found in the imaginary part. Whereas 
magnitudes of the imaginary parts in (I) and (II) are much larger than those in Fig. |6l they 
are suppressed in (III) to give moderate enhancement of the imaginary part at B r = 5000, 
as in the real part. These results indicate that stronger distortion of the Dirac sea induces 
stronger effects of the competing back reactions. As for the energy dependence, the stable 
branch extends up to infinite photon energy as in Fig. [6j The connecting point between 
the stable and unstable branches moves to a larger photon energy compared to the case in 
Fig. [6j indicating a larger upward shift of the threshold. The imaginary part has a broader 
profile in Fig. [7J We will investigate independence more systematically for the refractive 
index in Sec. IV Bl and the dielectric constant in Appendix [A] 



V. REFRACTIVE INDEX IN THE LLL APPROXIMATION 

Thus far, we have considered only the dielectric constants e. However, in some cases, it is 
more convenient to treat the refractive indices n, which can be immediately obtained from 
the dielectric constants through a simple relation n 2 = e. In this section, we will present 
results for the refractive index in the LLL approximation, and in particular we will focus on 
its dependences on the angle and magnetic field strength. 

When the dielectric constants have imaginary parts, e = e rca i + ze imag , we can similarly 
define real and imaginary parts of the refractive indices as n = n rea i + i n; mag , which are 
explicitly given with respect to the refractive index as, 

™rcal = -^=\/l e l + e roal, (19) 
1 



^irnag = -^=V|e| - Creal, ( 20 ) 



with the magnitude |e| being defined by |e| = ye 2 cal + ef mag . 

As we saw in the previous sections, only the parallel mode of the dielectric constant ey 
remains nontrivial in the LLL approximation, and the same is true for the refraction index: 
we have only n\\ as a nontrivial mode. Figured] shows the refractive index nn at B r = 500 as a 
function of (scaled) photon energy squared a) 2 that corresponds to the dielectric constant in 
Fig. El We find that the refractive index has a similar structure in photon energy dependence, 
and that it becomes complex when the dielectric constant takes the unstable branch. We do 
not repeat the explanation for the transition of the solutions from the naive prescription (I) 
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FIG. 8. Refractive index nu for 6 = 7r/4 and B r = 500: each line shows the corresponding 
quantities explained in the legend in Fig. [6j 



to the fully self-consistent treatment (III). All the explanations for the dielectric constant 
ey equally apply to the refractive index n\\. In this section, we rather discuss dependences 
of the self-consistent solutions on the propagation angle 9 and the magnetic field strength 
B T . In fact, the same analyses were in advance performed for the dielectric constant en as 
summarized in Appendix |A], and we have used those results to compute the refractive index 
through the relations f|T9l) and (1201) . 
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FIG. 9. Refractive index with varying propagation angle 8: The refractive index takes the same 
value, when photon propagation is oriented in the directions, ±9, ir ± 8. 

A. Dependence on propagation angle 




In the previous section and in Fig. [Sj, the results of ey and tin were shown only for the angle 
9 = 7r/4. However, there is a significant angle-dependence in the results, which we are going 
to present in this subsection. Upper and lower panels in Fig. |H] represent real and imaginary 
parts of the complex refractive index when the propagation angle varies from 9 = to tt/2. 
We have shown the results only in this region of 9 because the refractive index has the 
following symmetry: e\\(9) = e\\(—8) = e\\(7i ± 9). In the upper panel, solid lines are stable 
branches with different photon energies. A long-dashed line indicates a limiting value n\ [m of 
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the refractive index on the stable branch. It is given by num = y/e^ = | cos#| _1 , namely a 
square root of the limiting value of the dielectric constant ( [15]) . As seen in Fig. [8]for 9 = 7r/4, 
the refractive index approaches the limiting value nu m as the photon energy increases. At 
a certain photon energy above around u 2 ~ 1, the refractive index on the stable branch 
strongly depends on the angle. The real parts grow with increasing 9, indicating that, when 
a photon propagates almost perpendicularly to the magnetic field, the effects of vacuum 
birefringence appear more strongly and photon's phase velocity becomes significantly small. 

While we obtained the simple expression of the limiting value nu m within 1-loop accuracy 
of the vacuum polarization tensor, it is still an open question how higher-order diagrams 
contribute to the refractive index. Indeed, we notice that the limiting value rin m diverges 
in particular two cases when the photon propagates in the perpendicular directions 9 = 
±7r/2 with energy Co 2 — 1. These singular behaviors would suggest that further careful 
investigation is required when the modification of the refractive index is very large, e.g., the 
vacuum polarization would be suppressed by mutual Coulomb interaction between a fermion 
and antifermion pair appearing in the loop part which corresponds to photon exchanges in 
the multi-loop diagram. 

There are also unstable branches (short-dashed lines) when the angle is large. This is 
clearly seen in the lower panel where the imaginary parts are shown. The stable (unstable) 
branches appear in small (large) 9 region because the threshold condition r| = 1 depends on 
the angle. Also, the angle where unstable branch starts to appear depends on the photon 
energy, and decreases with increasing energies. Therefore, if the photon energy is large 
enough, it can decay even at small angles. Note, however, that the decay never occurs at 
9 = within constant magnetic field case, where the refractive index persists to be unity as 
in the ordinary vacuum. 



B. Dependence on the magnetic- field strength 

In this subsection, we present the dependence of the refractive index on the magnitude 
of the external magnetic field. The real part of the refractive index on stable and unstable 
branches, and the imaginary part on the unstable branch are respectively shown below. 
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FIG. 10. B r dependence of the real part of nu in the LLL approximation. Lines show the 
refractive index at 9 = tt/4 and various photon energies. Upper and lower panels show the stable 
and unstable branches, respectively. In the upper panel, results from the LLL approximation agree 
well with those from the numerical integration which are indicated with dots. 

1. Real part 

Upper panel of Fig. [10] shows a Independence of the refractive index nil on the stable 
branch. Since the stable branch is in the region below the threshold, the refractive index has 
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only a real part. Curves and dots indicate the results from the LLL approximation and the 
numerical integration, respectively, which agree well with each other in the strong-field limit 
B r > 10. Thus, we again confirm the validity of the LLL approximation when the magnetic 
field is strong enough B r > 10. We find that, at any photon energy, magnitude of the 
refractive index increases with increasing B r , and saturates in the limit of strong magnetic 
field. The limiting value is given by n Um = y/e liin as mentioned in the last subsection. While 
we discussed this limiting value n Um as reached in the limit of a large photon energy (see 
Eq. (I18p ). the plot shows the same limit is reached, at any photon energy, in the strong-field 
limit. Remind that the limiting value of the dielectric constant e,. is obtained in the limit 
in which the scalar function Xi becomes divergently large Xi ^ 1 as the photon momentum 
approaches the threshold ?" 2 (e||) — > 1. Note also that the scalar function is proportional 
to the field strength, xi <X- B T , showing a divergently large value also in the strong field 
limit (see also Fig. [4]). Therefore, even at a small photon energy, the dielectric constant 
approaches the limiting value e lim when the magnetic field is strong enough. 

In the limit of zero frequency Co 2 = 0, both r\ and r 2 are vanishing, and the scalar 
coefficient function Xi LL has a simple form (see Eq. OH])), 

X^ LL (0) = ^B v , (21) 

which indicates a monotonic increase with respect to B T . Thus, x\ hh becomes of the order 
x\ LL ~ 1 if the magnetic field is as strong as B T ~ 3n/a ~ 10 3 . Beyond this field strength, xi 
starts dominating unities in the denominator and numerator in Eq. (1T21) . Then, the dielectric 
constant and thus the refractive index would start approaching the limiting values e lim and 
n lim , respectively. The upper panel in Fig. ITUl shows sizable modification of the refractive 
index already around B r ~ 10 3 at zero and low energies {B r = 37r/a gives n» = 4/3 when 

e = tt/4). 

As the photon energy becomes large beyond u 2 = 1/(1 — cos 2 9) (a) 2 = 2 when 9 = tt/4 
in Fig. [TOi) . the refractive index is subject to significant modification even in a relatively 
weak regime, B r < 100. This is because the squared momentum approaches the threshold 
r|(e||) = 1, and thus Xi increases divergently (see Fig. [3]). 

Lower panel in Fig. [10] shows the real part of the refractive index on the unstable branch. 
In this plot, each line has an end point at a certain B r , because the unstable branch does 
not exist below a connection point between the stable and unstable branches, as shown in 
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FIG. 11. B T dependence of the imaginary part of nn on the unstable branch. Curves show those 
at 9 = 7r/4 and various photon energies. 



Fig. |8j Black long-dashed curve shows location of the connection point specified by a pair 
of photon energy and magnetic field strength. As found in a comparison between Fig. [6] 
and Fig. [7] the connection point shifts to a larger photon energy as B r increases. Therefore, 
for a given photon energy, unstable branch does not exist above a certain magnetic field 
strength. Although peaks apparently look sharper as the photon energy increases, this is 
due to a logarithmic scale for the horizontal axis. 



2. Imaginary part 

Figure [11] shows B T - dependence of the imaginary part of the refractive index on the 
unstable branch. As in the case of the real part shown in the lower panel in Fig. [10] there is 
a peak structure just below the end point of each curve, reflecting a prominent photon decay 
rate near the connection point. As we found in the comparison between Fig. [6] and Fig. the 
threshold obtained in the self-consistent solution (the connection point) shifts more distantly 
from the naive one u) 2 h = 1/ sin 2 9 to a large photon energy when the magnetic-field strength 
is larger. Consistently to what we found there, a stronger magnetic field gives rise to the 
peak for an energetic photon in Fig. [10] 
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FIG. 12. Trajectories in a complex ray-plane showing its magnitudes on the unstable branch: 
The propagation angles are commonly fixed at 9 = 7r/4, and the magnetic field strength on each 
trajectory is taken as B T = {50, 100, 500, 1000, 5000, 10000, 50000} from the innermost to the out- 
ermost. On each curve, the photon energy increases from the right to left edges, where the right 
edge corresponds to the connection point between the stable and unstable branches. 

3. Real vs imaginary part on the unstable branch 

Figure [12] shows a relation between real and imaginary parts of the complex refractive 
index on the unstable branch. The magnitude of the external magnetic field is fixed on each 
arc, and a radius increases as the magnetic field strength becomes larger. A gradation of 
colors transits from red to violet as the photon energy increases, of which scale is fixed so 
that the full range from red to violet appears on the innermost arc. Every arc has a left 
edge located on en = 1 + • i, because the refractive index on the unstable branch converges 
to unity in case of the strong-field limit within the LLL approximation, when the photon 
energy becomes large beyond the first threshold and stays small enough not to approach 
the second threshold. The opposite edge corresponds to the connection point between the 
stable and unstable branches. The imaginary part of the refractive index vanishes at this 
point, and the real part becomes large as the field strength becomes large. Magnitudes of 
the modification of the real and imaginary parts are comparable to each other at a fixed B T . 
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FIG. 13. Decay length of photon propagation in strong magnetic field: Photon decays within 
order of 0.1 - 10 picometer above the threshold, while magneter has magnetic field extending in a 
macroscopic scale. 

C. Decay length 

Lastly, let us show the decay length defined in Eq. (j3J). Figure [13] shows the decay 
length at a magnetic field strength B r = 500 and propagation angles, 6 = tt/6, 7r/3, tt/2. 
Photon energy at the decay threshold becomes small as the angle increases, because of an 
angle dependence of the threshold condition, r| = 1. While "magnetar" is thought to be 
accompanied by a strong magnetic fielcjl of order B r ~ 100 extending in a macroscopic scale, 
plotted lines show that photon decays within order of 0.1 - 10 picometer if photon energy 
is above the threshold. This indicates a drastic modification of gamma-ray and dilepton 
spectra emitted from magnetars. 



VI. SUMMARY AND PROSPECTS 



We investigated the vacuum birefringence phenomena in strong magnetic fields on the 
msis of the analytic representation of the vacuum polarization tensor obtained in paper I 
l|. In particular, we analyzed in detail the phenomena in the lowest Landau level (LLL) 
approximation where both the dielectric constant (en) and the refractive index (tin) of a 

7 Here, physical scale is provided by the critical magnetic field strength defined by electron mass. 
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propagating photon in the mode parallel to the external magnetic field deviate from unity 
and become complex. The LLL contribution corresponds to the first term of the polarization 
tensor which is represented as the infinite sum with respect to the Landau levels, and turned 
out to be a very good approximation when the magnetic field is strong enough B/B c > 1 
with B c being the critical magnetic field strength for the relevant fermions. 

In the LLL approximation, there is a threshold in the photon energy beyond which both 
en and nil aquire imaginary parts and thus the decay of a real photon into a pair of a fermion 
and an antifermion (in the LLLs) becomes possible. With the explicit analytic expression 
for the polarization tensor, we were able to thoroughly inspect properties of ey and n\\ such 
as the dependences on the propagating angle and the magnetic-field strength, and also how 
the imaginary parts appear. 

We found that the equations that define e (or equivalently n) with respect to the scalar 
coefficient functions in the polarization tensor are implicit functions of e (or n), and must be 
solved self-consistently in terms of e (or n). This self-consistent procedure is indispensable 
for accurate description of e and n when the deviation of e and n from unity is large, which 
is realized in the limit of strong fields or high-energy photons. Physically, this procedure is 
necessary for consistently taking into account the back reaction from the distorted Dirac sea 
in response to an incident photon field. Namely, an induced polarization and the photon 
decay bring in modifications of the external momentum, and thus dispersion, of the incident 
photon. 

The refractive index amounts to a large value in the strong field limit. For example, at 
the angle 9 = 7r/4, the refractive index n\\ keeps increasing with stronger magnetic fields, 
and approaches the limiting value nn m = \/2. This is comparable to the values which we 
encounter in ordinary life. To name a few, atmosphere of the earth (1 atm, °C) and 
water (20 °C) have refractive indices, n air = 1.000293 and n wa t er = 1.333, respectively, and 
"calcite" known as a representative birefringent material has refractive indices n Q = 1.6584 
and n e = 1.4864 for ordinary and extraordinary modes, respectively. At a magnetic field 
strength B r ~ (m^ / m e ) 2 ~ 10 4 " 5 which could be realized in the ultrarelativistic heavy-ion 
, the refractive index is close to the limiting value ny im = \/2, and thus can 
be larger than that of gas, and even comparable to those of liquid and solid. 

An accurate description of the complex refractive index across the lowest threshold will 
8 and m e are masses of a pion m w ~ 140 MeV and an electron m e ~ 0.5 MeV, respectively. 
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be directly applicable to the studies of intriguing phenomena in magneters whose magnetic 
field is estimated to be two orders larger than the critical magnetic field strength [la \v\ . 
In particular, descriptions of both the real and imaginary parts will work as robust building 
blocks necessary for investigating the interplay and competition among photon splitting, 



vacuum birefringence and photon decay 18|, [19 ] . 



Besides, as mentioned above, the ever strongest magnetic field will be created in the 
ultrarelativistic heavy-ion collisions in RHIC and LHC experiments. Since the magnitude 
of the magnetic field could be four orders larger than the critical magnetic field, radiations 
from the created matter in collision events would interact with the extremely strong magnetic 
fields, and bring out the information of the early-time dynamics 



20 



2l|| . However, photons 



created in heavy-ion collisions have, in general, large energies. Thus we need to extend the 
present work towards including important contributions among all the Landau levels shown 
in paper I. 

In application to laser physics 22j, it will be very important to include effects of a strong 
electric field as well as a magnetic field. For this sake, basic concepts and techniques exam- 
ined in this series of papers could be useful to find differences out of similarities between 
expressions of the vacuum polarization tensors in the presence of external electric and mag- 
netic fields. 
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Appendix A: Dependence of the dielectric constant on B T and 8 

In this supplementary section, we show dependences of the dielectric constant e\\ on the 
propagation angle 9 and the magnetic-field strength B T . By self-consistently solving Eq. (fT2"j) . 
we obtained plots shown in Figs. [M] - [TTJ, which are then mapped to the refractive index 
thought the relations (fT9|) and (|20|) to yield Figs. [9]-[T2j Parameters in these plots are taken 
as the same values as those in the corresponding plots of the refractive index in Fig. l9l-[T2l 
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FIG. 15. B r dependence of the real part of en on the stable (left) and unstable (right) branches. 
Photon's propagation angle is fixed at 9 = tt/4, and photon energies are varied as indicated in the 
panels. Dielectric constant on the stable branch, obtained from the numerical integration, is shown 
with dots in the left panel. 



The dependence on the photon's propagation angle is shown in Fig. HH The real part of 
the dielectric constant on the stable and unstable branches are shown in the left and right 
panels in Fig. [TJ)1 respectively, while the imaginary part on the unstable branch is shown 
in Fig. [T6j Magnitudes of the real and imaginary parts on the unstable branch are shown 
together in Fi.g[T7]as contours in the complex ey-plane. 
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Re[e„] 



FIG. 16. S r -dependence of the imaginary part FIG. 17. Trajectories on a complex ey-plane. 
of e 1 1 on the unstable branch. Parameters are Parameters are the same as in Fig. [121 
the same as in the right panel of Fig. [15j 

Appendix B: Remark on other Landau levels 



We remark on structures near the thresholds in other Landau levels which have the same 
property as that of the lowest Landau level. In the infinite series expression of the scalar 
coefficient functions x« with all the Landau levels (see Sec. IV in paper I [l|), we find a 
relation between the coefficient functions xo — ~X2, if we pick up an individual contribution 
specified by i = and 0). It immediately follows from Eq. ([T|) that the dielectric 
constant of the perpendicular mode is then given by 



1 + Xo 



1 + Xo cos2 # 



(Bl) 



The right-hand side of the above expression has the same form as Eq. ( 112]) . except that the 
scalar function xi is replaced by xo- Because Xo diverges at these thresholds as Xi does at 
the lowest threshold, the dielectric constant has a liming value, 



lim e_i_ (rf<) 



1 



cos 2 9 



(B2) 



when the photon momentum approaches the thresholds as rjj — > (1 + y/T+ 2nB r ) 2 /A = s ^. 
The right-hand side of Eq. (IB 2 1) is exactly the same as that of Eq. ( fl5l) . Therefore, the 
dielectric constant e± in the perpendicular mode would, in the region around the thresholds 
at £ = and 0), have the same structure in the photon energy dependence as that of 
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the longitudinal mode en around the lowest threshold. 
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